Electromagnetic Leptogenesis 
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We present a new leptogenesis scenario, where the lepton asymmetry is generated by CP violating 
decays of heavy electroweak singlet neutrinos via electromagnetic dipole moment couplings to the 
ordinary light neutrinos. Akin to the usual scenario where the decays are mediated through Yukawa 
interactions, we have shown, by explicit calculations, that the desired asymmetry can be produced 
through the interference of the corresponding tree-level and one-loop decay amplitudes involving 
the effective dipole moment operators. We also find that the relationship of the leptogenesis scale 
to the light neutrino masses is similar to that for the standard Yukawa-mediated mechanism. 
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I. INTRODUCTION 

Baryogenesis via thermal leptogenesis [l[ provides an 
elegant explanation of the cosmic baryon asymmetry [2| • 
The conventional setup involves minimally extending 
the Standard Model (SM) by adding three heavy right- 
handed (RH) Majorana neutrinos, which are electroweak 
singlets, and allowing them to interact with ordinary left- 
handed (LH) lepton doublets via complex Yukawa cou- 
plings. As a result, when the heavy neutrinos decay out- 
of-equilibrium in a CP violating way, a lepton asymme- 
try is generated in the early Universe. This asymmetry 
is then partially converted by sphaleron processes to the 
baryon asymmetry we detect today. 

An attractive aspect of this scenario is that it naturally 
allows the implementation of the see-saw mechanism [3[ 
which gives the light neutrinos a tiny but nonzero mass. 1 
Consequently, a strong link between some neutrino prop- 
erties and successful asymmetry generation can be es- 
tablished. For instance, to simultaneously obtain a suffi- 
ciently large lepton asymmetry and the correct light neu- 
trino masses, the heavy neutrino masses must be larger 
than 10 9 GeV in most viable leptogenesis scenarios [yj. 

In this paper, we will consider a scenario where leptoge- 
nesis is mediated not by the standard Yukawa couplings, 
but instead by electromagnetic dipole moment couplings. 
One motivation for introducing such couplings is to ex- 
plore whether new sources of CP violation can lead to a 
significant lepton asymmetry. A natural question we may 
ask is whether the introduction of CP violating dipole 
moment couplings will allow leptogenesis to occur at a 
lower scale, closer to experimentally accessible energies. 
However, we find this not to be possible due to the con- 
nection between the dipole moment couplings and the 
neutrino mass. 

The general form of a dipole moment coupling of the 
light neutrinos, v, to the heavy neutrinos, N, is given by 
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T7(/i + id r yc,)cr a PNF a p, where n and d are the magnetic 
and electric transition moments, respectively. These 
dimension- 5 effective operators may be assumed to be 
generated by some new physics beyond the electroweak 
scale. While we do not speculate on the nature of this 
new physics, the inclusion of these operators permits a 
new leptogenesis mechanism. Radiative decays of the 
heavy neutrinos, N — > v + 7, can now produce the re- 
quired lepton asymmetry in the early Universe, provided 
that the complex electromagnetic dipole moment cou- 
plings violate CP. 

Below, we outline the relevant properties of the elec- 
tromagnetic dipole moment (EMDM) couplings, and dis- 
cuss the necessary requirements for a decay process to 
manifestly violate CP. We shall then explicitly calculate 
the decay rates and CP asymmetry for electromagnetic 
leptogenesis, and compare with the standard Yukawa- 
mediated leptogenesis scenario. 



II. ELECTROMAGNETIC COUPLING 
BETWEEN LIGHT AND HEAVY NEUTRINOS 



We extend the particle content of the minimal SM by 
adding heavy RH neutrinos, Nr, which are assumed to 
have large Majorana masses. Since we are interested in 
leptogenesis energy scales above the electroweak phase 
transition, we will take the usual light neutrinos to be 
massless LH states. The most general electromagnetic 
dipole moment coupling of the heavy RH neutrinos to 
the light LH neutrinos is then given by 



T-Cem = gijVLi a a/3 N Rj F a p + h.c. 



(1) 



1 Note that similar results can be instead achieved with the addi- 
tion of heavy Higgs triplets 4] , or with other mechanisms [j| . 



where g is a complex (dimensionful) matrix, and 
are flavor indices. Note that there is only one distinct 
electromagnetic dipole moment coupling when expressed 
in term of LH and RH chiral fields (rather than dis- 
tinct magnetic and electric dipole moment terms) since 
75-Pl,_r oc P L ,fl- 
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III. CP VIOLATION IN DECAYS 



IV. A TOY MODEL 



If the dipole coupling of Eq.fTT]) had less-constrained 
chiral structure, then for a given Vi and Nj there could be 
independent magnetic and electric transition moments, 
Hij and dij . It might be thought that in this situation, a 
tree-level interference between the amplitude induced by 
Hij and that induced by could lead to a CP-violating 
difference between the rates for Nj — > v% + 7 and its CP 
conjugate. However, there can never be a difference be- 
tween the rates for CP-conjugate decay modes until one 
goes beyond first order in the underlying Hamiltonian. 
This fact is well known, but it is interesting to see that 
it can be proved very simply by using CPT invariance. 
Consider, in the rest frame of the parent particle Q, the 
decay Q — > a% + a 2 + • • ■ . If CPT invariance holds, the 
amplitude for this decay obeys the constraint 

\(a 1 (p 1: X 1 )a 2 (p 2 ,X 2 )---\T\Q(m))\ 2 

= |(ai(pi,-A 1 )a 2 (p 2 ,-A 2 ).-.|Tt|Q(-m))| 2 . (2) 

Here, pi and Aj are, respectively, the momentum and 
hclicity of daughter particle Ui, fh is the z-axis projection 
of the spin of Q, and T is the transition operator for the 
decay. If S is the ^-matrix operator, T = i(S — I)- To 
first order in the Hamiltonian Ti for the system, T = Ti, 
so that, to this order, T* = T. From the latter relation 
and Eq.Q, it follows that, after summing over the final 
helicities and integrating over the outgoing momenta, 



r [Q — > ai + 0,2 -I ] = T [Q — > a\ + a 2 + 



(3) 



This equality must hold to first order in Ti regardless of 
whether Ti contains numerous terms and CP- violating 
coupling constants. 

In the special case of a two-body decay, Q — > a\ + a 2 , 
we have pi = —p 2 = p. For this case, let us rotate the 
system of particles on the right-hand side of Eq.@ by 
180° about the axis perpendicular to the z-axis and to 
p. Eq.([2]) then states that, to first order in Ti (so that 

r = rt), 

|(ai(p,Ai)a 2 (-f,A 2 )|r|Q(m))| 2 

= |<ai(-£-Ai)a 2 (a-A 3 )|T|3(a))| 2 . (4) 

The processes whose amplitudes appear on the two sides 
of this constraint are the CP- mirror images of each other. 
Thus, in two-body decays, to first order in Ti, the rates 
for CP-mirror-image decay processes must be equal even 
before one sums over final helicities and integrates over 
outgoing momenta. 

We conclude that CP-violating rate differences be- 
tween CP-conjugate electromagnetic decays of heavy 
neutrinos can only appear once amplitudes involving 
loops are included. 



In this section we illustrate, by means of a toy model, 
the viability of generating a lepton asymmetry through 
EMDM interactions between ordinary light neutrinos and 
the postulated heavy Majorana neutrinos in the early 
universe. In order to illustrate the physics as transpar- 
ently as possible, we begin by considering a simplistic 
model which is not invariant under the SM gauge symme- 
try, SU(2)l x U(1)y, but is instead invariant only under 
the electromagnetic symmetry U(1)q. We will general- 
ize to a realistic model in which invariance under the SM 
gauge group is enforced in section IVl below. 

We assume the EMDM couplings are generated by new 
physics at an energy scale A > M, where M denotes a 
heavy Majorana neutrino mass. We work with an effec- 
tive theory that is valid below the scale A, obtained after 
integrating out all new heavy degrees of freedom. The 
lowest dimension EMDM operator of interest in such a 
scenario is given by: 



r 5B 



= —e-^Xjk; 



V Lj cr a Pp R N k F a p + h.c. , (5) 



A 



where j = e,fi,T and k = 1,2,3. The electromagnetic 
field strength tensor is F a p = d a Ap — dpA a , with A a 
being the photon field. The light neutrinos are denoted 
by Vj, while Nk is the heavy neutrino field in the mass 
eigenbasis which satisfies the Majorana condition: N k = 
e l(f>k N^, for some arbitrary phase cpk- We have defined 
A as a dimensionless 3x3 matrix of complex coupling 
constants, while A is the cut off scale of our effective 
theory. For convenience, we have factored out e~ l(fk ^ 2 
for each A'fc in Eq.© so that the Majorana phases will 
not appear explicitly in any of our final expressions. 

To ascertain whether leptogenesis is possible, the key 
quantity of interest is the CP asymmetry in the decays 
ofN k : 



-(5) 
'-k,j 



(N k ^v 3 7) 



r 



' r (JV fe ->F7) 



(6) 



where r^^;^) = J2j ^(N k ^Uj 7)- The lowest order con- 
tribution to the decay rate, shown in Fig. [2 is given by 



= r 



(AU) fcfc M, 



4tt 



— • (7) 

A 2 [ 1 



(5) 

The leading contribution to the CP asymmetry, e k p 
comes from the interference of the tree-level process of 
Fig. [1] with the 1-loop diagrams with on-shell intermedi- 
ate states depicted in Fig. [2] As with standard (Yukawa 
interaction) leptogenesis, the CP asymmetry receives two 
contributions: the self-energy and vertex correction. We 
have calculated these explicitly and obtain 



r (5) 

"self-fejj 



-(M fc /A) 2 v 
(A t A)fe m ^~ 



Im 



(AU) r 



(8) 
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FIG. 1: The tree- level diagram for the decay of N k via the 
EMDM interaction of Eq. ©. 




(a) 



(b) 



FIG. 2: (a) Vertex and (b) self-energy diagrams which con- 
tribute to the CP asymmetry of Nk decay via the interaction 
of Eq.([5]). Note that since weak isospin is violated in this 
model, both v n and V n are allowed in the loop of (a). 



for the self-energy contribution and 
-(M fc /A) 2 



J5) 



27r(AU) fefe 



£ Im [\* k \ jm ^X) km f(z)] 



(9) 



with f(z) = \[z 



1 + 2z 1 - (z + 1) In 



"z+r 




z 





for the vertex piece, where 2 = M^/M?. Note that in 
Eq.© and Eq.©, we have not yet summed over the final 
lepton flavor j . 

The expressions given in Eq.© and Eq.© are very 
akin to those in standard leptogenesis [7[. The first 
and second terms of Eq.© correspond to the interfer- 
ence terms involving 77„ and v n as the intermediate state 
in Fig. b) respectively. It should be emphasized that 
upon summing over j (ie. ignoring flavor effects [1]), the 
second term vanishes. Furthermore, explicit calculations 
have shown that in Fig.[2ja), the contribution from the v n 
intermediate state actually evaluates to zero, and hence 
there is no term proportional to (A^A) m fc in Eq.©. 

The total CP asymmetry from the decays of N k s into 
light neutrinos and photons in this model is simply given 

b y e S = E 2f-fc,j + ffcLj- This asymmetry would be 
nonzero as long as there are phases in the coupling ma- 
trix A which cannot be removed by redefinitions of the 
neutrino fields. As A is an arbitrary complex matrix, it is 
not hard to see that one cannot eliminate all the relevant 
phases to render both Eq.© and Eq.© zero. Hence, 
this type of EMDM interaction between light and heavy 
neutrinos will in general generate a lepton asymmetry in 
the early Universe. Before discussing the magnitude of 
this asymmetry, we will first generalize this scenario to 
one in which the EMDM couplings respect the SM gauge 
symmetries. 



V. A MORE REALISTIC EXTENSION 

While the simplistic model in section IIVI can demon- 
strate the viability of lepton generation through EMDM 
operators, it is nonetheless unrealistic as it is incompat- 
ible with the SM. We now overcome this by considering 
only EMDM type operators that respect the SM gauge 
group. Again, we construct an effective theory by taking 
the usual minimally extended SM Lagrangian with three 
generations of heavy Majorana neutrinos, and augment- 
ing it with EMDM operators. The most economical of 
such operators involving only (the minimally extended) 
SM fields are of dimension six Q, and the interaction 
Lagrangian of interest is 



C 



EM 



1 



a/3 



P R N k + h.c. 



(10) 



where the r, are the SU(2)l generators, tj 
is the lepton doublet, and cj> = (</> *, — _ ) T is the SM 
Higgs doublet. The field strength tensors of U(l)y and 
SU(2) L are given by B af3 = d a Bp - dpB a and = 
d a Wp - dpWa - ge im nW™Wp , respectively, where g' and 
g are the corresponding coupling constants. Again, A is 
the high energy cut off of our effective theory, while the 
matrices of dimensionless coupling constants, A' and A', 
are in general complex. 

The higher dimension (non-renormalizable) operators 
of Eq. (fT0|) are assumed to be generated at the energy 
scale A, beyond the electroweak scale. The presence of 
these operators would imply the existence of some new 
physics at a high energy. After SU{2)l®U{1)y breaking, 
these operators will give rise to electromagnetic dipole 
transition moments of N and v. However, for the pur- 
poses of leptogenesis, we are interested here in the regime 
above the electroweak symmetry breaking scale. 

The decay of TV will now produce 3-body final states, 
namely N k — > £j(/)W^ and Nk — > £j(f>B a , as shown in 
Fig-E^a). Likewise, the self-energy and vertex corrections 
now become two-loop processes, an example of which is 
shown in Fig. [3Jb). As before, a lepton asymmetry will 
be generated through the interference of the tree-level 
amplitude with the on-shell part of the vertex and self- 
energy amplitudes. 

To demonstrate that the EMDM interactions in this 
model can indeed generate a lepton asymmetry, we 
have explicitly calculated the tree-level decay rate, 
^ {N k ^t.j<t>B a )i and self-energy contribution to the CP 

asymmetry, £^} f _ k j , for the couplings which involve the 
hypercharge boson B (see Fig. [3]). The SU{2) gauge bo- 
son will make a contribution to the decay rate and CP 
asymmetry of similar magnitude. For simplicity, we have 
not explicitly calculated the vertex corrections, an exam- 
ple of which is shown in Fig. 21 The vertex diagrams 
can again have on-shell intermediate states and, barring 
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(a) (b) 

FIG. 3: (a) The tree-level diagram for the the 3-body decay: 
N k -> £j(j>B a induced by the first term in Eq.JlO}. (b) The 
corresponding self-energy diagram. 




FIG. 4: An example of a vertex correction to the 3-body decay 
process. 



JV 




v 



FIG. 5: Contribution to the neutrino Dirac mass, ran, in- 
duced by the electromagnetic dipole moment operator. 
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FIG. 6: Contribution to the light neutrino mass, , induced 
by the electromagnetic dipole moment operator. 



accidental cancellations, would contribute to the lepton 
asymmetry generated. 

For the diagrams shown in Fig. [3j our explicit calcu- 
lations of r and e lead to similar forms to those found 
earlier: 2 



and 



-(6) I 
-sclf-fc ,j I 



~ (H) r < w *-Vr) > 



( M k 



\8ttA 



(5) 

sclf-fc, j I 1 



(11) 



(12) 



where we must replace A — > A' in the RHS of Eq. (TTTj) and 
Eq. (fl~2"]) . It is thus possible to generate a non-zero CP 
asymmetry via the EMDM type interactions of Eq. tfTO"]) . 
provided A' contains complex phases. 



VI. NEUTRINO MASS 

We now address the connection between the dipole mo- 
ment operators and the neutrino mass. The new physics 
that gives rise to the effective operators in Eq. (JTUJ) might 
also be expected to give rise to neutrino mass terms. It is 
well known that via a careful choice for the new physics, 
one can obtain large neutrino dipole moments without 
correspondingly large mass terms [ill, E2, EH E3, EH- 
However, even if neutrino masses are absent or sup- 
pressed at lowest order, radiative corrections involving 
the dipole moment operators generically induce neutrino 
mass terms [IEEE!- 



2 It should be noted that these 3-body decay processes are inher- 
ently different from those discussed in [Hj for these are 2-loop 
rather than 1-loop diagrams. 



As in Ref. 0], the operators in Eq. (fTTj)) will lead to a 
Dirac mass term for the neutrinos via the diagram shown 
in Fig. [5l There is no model-independent way of calcu- 
lating the exact size of this neutrino mass contribution, 
since the exact relationship between the dipole moments 
and the mass requires a UV completion of the theory 
(i.e. it depends on the nature of the physics at scale A). 
However, we may estimate the size of the contribution to 
the neutrino mass using naive dimensional analysis. The 
neutrino Dirac mass arising from Fig. [5] is thus estimated 
to be 



A' 



m D 



A 2 16tt 2 



WA 2 



(13) 



where the A 2 in the numerator arises from the cut off of 
the loop integral, and g' is the gauge coupling constant. 
This Dirac mass term will lead to a contribution to the 
light neutrino mass via the see-saw mechanism of 



m D M mo 



fJ 



16tt 5 



(14) 



However, there will also be a direct contribution to Majo- 
rana mass of the light neutrinos via the diagram in Fig. [51 
which we estimate as 



m 



B - \' T MX 



1 



A 2 16tt 2 ' 



(15) 



In the standard Yukawa-mediated leptogenesis sce- 
nario, the light neutrino masses are linked with the lepto- 
genesis parameters, since the Yukawa coupling constants 
that control the decay rate of the N also appear in the 
Dirac mass terms and thus (via the seesaw mechanism) 
in the v masses. In order for these Yukawa coupling con- 
stants to give rise to the correct values for both e and m„, 
the heavy neutrinos are required to have masses M > 10 
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GeV [6| 3 . The presence of the neutrino mass term aris- 
ing from the diagrams of Fig. [5] and Fig. [5] signifies that 
a similar connection will hold for the couplings A' that 
control electromagnetic leptogenesis. Hence, we will also 
need a relatively high N mass scale in order to produce 
a sufficiently large asymmetry together with sufficiently 
small v masses. 



Yukawa 


Electromagnetic 


Ti = IT- f h ) " M i 


r. 


1 Im(/ l t /i)? m Mi 


1 Im(A' T A')i m Mi / Mi 2 > * 


£ ~ 8tt (h*h)u M m 


6 ~ 2tt (A'tA')n M m V87rA 2 y 




d A' T MA' l2 1 
A* W"l8^ 



TABLE I: Comparison of key quantities in standard and elec- 
tromagnetic leptogenesis, where h denotes Yukawa coupling 
constants. We have assumed there is at least a mild hierarchy 
in the masses of the heavy neutrinos, such that the asymme- 
try is predominantly generated from the decay of the lightest 
state, Ni. 



VII. DISCUSSION 

If we assume there is at least a mild hierarchy in the 
masses of the the heavy neutrinos, the asymmetry will 
be predominantly generated from the decay of the light- 
est state, N\. In Table HI we compare the decay rate, 
CP asymmetry, and light neutrino masses for electro- 
magnetic leptogenesis with the corresponding expressions 
for the standard (Yukawa) scenario. The coupling con- 
stants and the N\ mass enter into the expressions for Ti, 
e, and m„ in essentially the same way for the two sce- 
narios. Modulo the additional suppression factors in the 
RHS of the table, we see that the region of viable pa- 
rameter space for electromagnetic leptogenesis must be 
similar to that for the standard Yukawa mechanism. It 
is also clear that although we require A > M in order 
for our effective operator approach to be valid, we do not 
want A 3> M, as it would suppress both T and e. We 



This condition may be alleviated in degenerate leptogenesis sce- 
narios, in which the CP asymmetry can be enhanced and thus 
the N mass lowered through a resonance in the self-energy con- 
tribution when Mfc ~ M m list . 



therefore require a moderate hierarchy between M and A 
in order to obtain an appropriately sized CP asymmetry. 

In general, both the Yukawa and EMDM interactions 
will contribute to the decay rate, CP asymmetry, and 
neutrino mass. Depending on the relative size of h and 
A', either one mechanism will dominate or there will be an 
interplay between the two. In what follows, we suppose 
that the Yukawa couplings are negligible. For simplicity, 
we also ignore the flavour structure of the matrices A', 
m v , and M, and make the crude assumption that all 
elements are of similar magnitude. If we then take 

A ~ 10M 2 ~ 20M X , and A' - 35, (16) 

we obtain an asymmetry of e ~ I0 -6 , while smaller val- 
ues of A' would lead to a correspondingly smaller asym- 
metry according to e oc (A') 2 . We define the decay pa- 
rameter as 

K = H\ T = Ml , (17) 

where H is the Hubble expansion rate. Note that K 
controls whether the N\ decays are in equilibrium, and 
is also a measure of washout effects via inverse decays. If 
we now take 

Mi - 5 x 10 12 GeV (18) 

together with Eq.([16|l. we obtain K — 0.3. Since K < 1 
(K 1) corresponds to the weak (strong) washout 
regime, this parameter choice should lead to moderate 
washout. Moreover, for these parameters the light neu- 
trino mass terms become ~ 0.04 eV and ~ 
0.1 eV, such that the neutrino mass is dominated by the 
contribution from the diagram in Fig. [5] For a larger 
Mi /A hierarchy, would be suppressed with respect 
to , but this is undesirable as the asymmetry e would 
also be suppressed by 4 powers of Mi/ K. 

Finally, we note that effective dipole moment interac- 
tions of two light neutrino states are induced by our La- 
grangian. The largest contributions arise from two-loop 
diagrams for which we estimate 

~ x - 5 x io_i9 ms ' (i9) 

where the second approximate equality assumes the pa- 
rameter values specified above. These induced dipole mo- 
ments are thus well below current experimental upper 
limits [19( which are of order 10" 11 /j,b- 

VIII. CONCLUSIONS 

In summary, we have presented a new leptogenesis 
mechanism in which electromagnetic dipole moment cou- 
plings induce CP violating decays of heavy RH neutri- 
nos. Via explicit calculation of the decay rates, we have 
demonstrated that a sufficient asymmetry can be pro- 
duced through such decays to make this scenario viable. 
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However, since the electromagnetic dipole moment oper- 
ators induce neutrino mass terms, the leptogenesis pa- 
rameters are constrained by the masses of the usual light 
neutrinos. For this reason, we find that leptogenesis must 
take place at a high mass scale, comparable to that for 
the standard scenario. 

Washout effects, which reduce the asymmetry, will dif- 
fer from the standard scenario. For example, the inverse 
decay is now a 3 — > 1 process, and there will be a differ- 
ent set of L violating scattering processes. We leave the 
detailed study of these washout effects to future work. 
Finally, it is possible that the interference of the tree- 
level and vertex amplitudes, which we have not explic- 
itly calculated, may enhance the asymmetry. As with the 
standard mechanism, a near degeneracy in the masses of 



the heavy neutrinos would also enhance the self-energy 
contribution to the asymmetry [3, HH . 
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